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Abstract — In recent work, Benjamin Scliumaclier and 
Miciiael D. Westmoreland investigate a version of quantum 
mechanics which they call modal quantum theory but which we 
prefer to call discrete quantum theory. This theory is obtained 
by instantiating the mathematical framework of Hilbert spaces 
with a finite field instead of the field of complex numbers. This 
instantiation collapses much the structure of actual quantum 
mechanics but retains several of its distinguishing character- 
istics including the notions of superposition, interference, and 
entanglement. Furthermore, discrete quantum theory excludes 
local hidden variable models, has a no-cloning theorem, and can 
express natural counterparts of quantum information protocols 
such as superdense coding and teleportation. 

Our first result is to distill a model of discrete quantum 
computing from this quantum theory. The model is expressed 
using a monadic metalanguage built on top of a universal 
reversible language for finite computations, and hence is di- 
rectly implementable in a language like Haskell. In addition to 
superpositions and invertible linear maps, the model includes 
conventional programming constructs including pairs, sums, 
higher-order functions, and recursion. Our second result is 
to relate this programming model to relational programming, 
e.g., a pure version of Prolog over finite relations. Surprisingly 
discrete quantum computing is identical to conventional logic 
programming except for a small twist that is responsible for 
all the "quantum-ness." The twist occurs when merging sets 
of answers computed by several alternatives: the answers are 
combined using an exclusive version of logical disjunction. In 
other words, the two branches of a choice junction exhibit an 
interference effect: an answer is produced from the junction if it 
occurs in one or the other branch but not both. 

I. The Result, Informally 
Consider this Prolog program: 

r ( false , false ) . 
r ( false , true) . 
r (true , false) . 

q(x) ;- r(false,X) . 
q(x) :- r(true,x) . 

The program starts with three facts about a relation r and 
then defines a rule q such that q ( X ) is true if either of the two 



clauses can be satisfied. Executing the query q ( X ) 
three answers: 



returns 



X = false ; 
X = true ; 
X = false. 



Now consider the same example expressed in the Discrete 
Quantum Theory over the field of booleans (DQT2) recently 
developed by Schumacher and Westmoreland |[l]: 



r|0) = 
r|l) ^ 



|0) 
|0) 



|1) 



q ^ r|0)+r|l) 

The relation r which could relate false to either false 
or true is expressed as a quantum gate r that maps |0) to 
a superposition of |0) and The rule q corresponds to 
a vector q produced by taking the superposition of the two 
possible alternatives. 

In contrast to the Prolog program, measuring the vector q in 
the standard basis {|0), |1)} is guaranteed to always return |1): 
the answer |0) could never be produced! To understand why, 
we calculate the result of q as follows: 

g = r|0)+r|l) = (|0) + |l)) + |0) = |l) 

In the last step the two intermediate answers |0} interfere 
destructively with each other and are canceled. The details 
are explained in Section |IV] 

This simple example captures the essence of our result 
which can be informally stated as follows. The model of com- 
putation inherent in DQT2, DQC2, is a relational program- 
ming model in which any value that appears an even number of 
times in disjunctions disappears. Section IVlII-CI shows a more 
significant example that implements the superdense coding 
quantum protocol in both DQC2 and Prolog. Removing all 
answers that appear an even number of times in the Prolog 
execution is indeed consistent with the quantum protocol. 

II. The Result, Formally 
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The pioneering research programmes of Abramsky and 
Coecke et al. 12, (|3l and of SeHnger 14j, Q have established 



that quantum computing (QC) can be modeled using dag- 
ger compact closed categories. The traditional mathematical 
model of QC — the category FDHilb of finite dimensional 
Hilbert spaces and linear maps — is a prime example of 
such categories, and so is the category Rely of sets and 
relations. Our result, formally stated, is that the computational 
structures inherent in DQC2 can be modeled in a non-standard 
category Reltti of sets and relations which is isomorphic to the 
category FDVec2 of finite dimensional vector spaces over the 
field of booleans and linear maps. This category is dagger 
compact closed and hence possesses all the computational 
structures necessary for QC. Our category differs from the 
standard category of sets and relations in one aspect; it uses the 
exclusive union of sets everywhere the standard union would 
be used. Diagrammatically we have: 

QC FDHilb 



DQC2 Rela/FDVec2 



Prolog Rel u 

III. Significance and Background 

Hilbert Spaces and Models of Quantum Computing: The 
traditional mathematical formulation of Quantum Mechanics 
(QM) is founded on Hilbert spaces. Although there are several 
other more abstract mathematical formulations of QM, the 
Hilbert space formalism remains the most accepted and the 
most widely used such formalism fE\. A Hilbert space is 
defined to be a real or complex inner product space that is 
a complete metric space with respect to the distance function 
induced by the inner product. A real or complex inner product 
space is a vector space over the field of real or complex 
numbers on which there is an inner product. In traditional 
QM, the underlying field is the field of complex numbers 
which serves as the space from which probability amplitudes 
for quantum events are drawn. 

For the purposes of QC, the mathematical formalism of 
QM is typically restricted to finite dimensional Hilbert spaces. 
This restriction removes one of the infinities in the formalism 
but retains another — the underlying infinite field of complex 
numbers. The resulting model is clearly still a mathematical 
idealization as it allows infinitesimally fine distinctions among 
quantum states that might differ by vanishingly small proba- 
bility amplitudes. In contrast, any particular realization of a 
quantum algorithm can only assume finite and discrete levels 
of representation of such probability amplitudes. Although it is 
customary for (classical or quantum) models of computations 
to include infinite structures of various kinds, it is important to 
understand how the idealized models emerge as the limits of 
finite approximations. A fundamental question that therefore 
motivates our research is whether it is possible to replace the 
infinite field of complex numbers by successively larger and 



larger finite fields to reach full QC in the limit. More impor- 
tantly, what is the computational power of these intermediate 
models, DQC„, as n grows larger? 

To tackle these questions, we do not work directly with the 
Hilbert space formalism, however. Indeed, the mathematical 
formulation of QC based on Hilbert spaces obscures many 
traditional computational structures that are inherent in the 
physical theory. Early on, Mu and Bird |7| showed that the 
simple model of QC based on vectors and linear operators is 
a monad which is the standard way to model computational 
effects |8|. In a subsequent development, the third author and 
collaborators L9J, LIOJ, [llj, \\2\, |13| have established that 
the more general quantum model based on density matrices 
and superoperators is an instance of the mathematical con- 
cept of arrows lfT4l which is a generalization of monads. 
Interestingly dagger closed compact categories identified by 
Abramsky et. al. also serve as models of Unear logic flS'l and 
various computational effects [16|, 1 17l . ITSll . Hence, based on 
these connections, our technical contributions are expressed 
using traditional computational structures and constructions: 
monads, arrows, and category theory. 

IV. Discrete Quantum Theory 

In their recent work, Schumacher and Westmoreland |1 1 ar- 
gue that much of the structure of traditional QM is maintained 
in the presence of finite fields. In particular, they establish that 
the quantum theory based on the finite field of booleans retains 
the following characteristics of QM: 

> the notions of superposition, interference, entanglement, 
and mixed states of quantum systems; 

• the time evolution of quantum systems using invertible 
linear operators; 

• the complementarity of incompatible observables; 

• the exclusion of local hidden variable theories and the 
impossibility of cloning quantum states; and 

• the presence of natural counterparts of quantum informa- 
tion protocols such as superdense coding and teleporta- 
tion. 

Fields: A field is an algebraic structure with notions of 
addition and multiplication that satisfy the usual axioms. The 
rationals, reals, complex numbers, and quaternions form fields 
that are infinite. There are also finite fields that satisfy the same 
set of axioms. Finite fields are necessarily "cyclic." 

We fix a field B consisting of two scalars {f, t}. The 
elements F and T are associated with the probabilities of 
quantum events, with F interpreted as definitely no and T 
interpreted as possibly yes^ The field B comes with an 
addition operation V (which in this case must be exclusive-or) 
and a multiplication operation A (which in this case must be 
conjunction). In particular we have: 

'Everything works if we switch the interpretation with F interpreted as 
possibly no and T as definitely yes [19] . 



F V F 

F y T 

T y F 

T V T 



F A F 
FAT 
T A F 
TAT 



The definitions are intuitively consistent with the interpre- 
tation of scalars as probabiHties for quantum events except 
that it appears strange to have TyT be defined as F, i.e., to 
have a twice-possible event become impossible. This results 
from the cyclic property intrinsic to finite fields requiring 
the existence of an inverse to y. This inverse must be T 
itself which means that T essentially plays both the roles of 
"possible with phase 1" and "possible with phase -1" and that 
the two occurrences cancel each other in a superposition. 

Vector Spaces: In the vector space over the field B, 
if vectors are represented as functions k that map basis 
elements v to scalars field values, then vector addition and 
scalar multiplication are defined by lifting the field operations 
to vectors as follows: 



Kl + K2 

b * K 



S.t. = Kl{v) y K2{v) 

s.t. k'{v) ~ b a k{v) 



Scalar multiplication is uninteresting as it either leaves the 
vector unchanged or returns the zero vector (denoted •) which 
maps everything to the scalar zero. Vector addition however is 
the key to introducing superpositions and interference effects. 

Consider the simple case of a 1-qubit system with bases |0) 
and The Hilbert space framework allows us to construct 
an infinite number of states for the qubit all of the form 
a|0) + /3|1) with a and /? elements of the underlying field 
of complex numbers and with the side condition that the 
state is not the zero vector and that it has length 1, i.e., 
that |ap + |/3p = 1. Moving to a finite field immediately 
hmits the set of possible states as the coefficients a and (3 
are now drawn from a finite set. In other words, in the 
field B, there are exactly three valid states for the qubit: 
f|0) + T|1) (which is equivalent to |1)), t|0) + f|1) (which 
is equivalent to |0)), and t|0) + t|1} (which we write as 
1+)). The fourth possibility is the zero vector which is not 
an allowed quantum state (discussed below). In a larger field 
with three scalars, there would be eight possible states for the 
qubit which intuitively suggests that one must "pay" for the 
amount of desired superpositions: the larger the finite field, the 
more states are present with the full Bloch sphere seemingly 
appearing at the "limit" n -> oo. 

Interestingly, we can easily check that the three possible 
vectors for a 1-qubit state are linearly dependent with any 
pair of vectors expressing the third as a superposition: 

|0) + |1) = 1+) 
|0) + 1+) = |1) 

|i) + 1+) - |0) 

In other words, other than the standard basis consisting of |0) 
and there are just two other possible bases for this vector 
space, {|1),|+)} and {|+), |0)}. 



The example also shows that the cyclic structure of the field 
extends to the vector space. 

Inner Products: A Hilbert space comes equipped with an 
inner product {vi \ V2) which is an operation that associates 
each pair of vectors with a complex number scalar value that 
quantifies the "closeness" of the two vectors. The inner product 
induces a norm \/ {v \ v) that can be thought of as the length 
of vector w. In a finite field, we can still define an operation 
{vi I V2) which, following our interpretation of the scalars, 
would need to return F if the vectors are definitely not the 
same and T if the vectors are possibly the same. This operation 
however does not yield an inner product, as the definition of 
inner products requires that the field has characteristic 0, i.e., 
that the repeated addition of T to itself never reaches F. As 
the above paragraph shows this is not the case for B (nor for 
any finite field for that matter) as the sum of positive elements 
must eventually "wrap around." In other words, if we choose 
to instantiate the mathematical framework of Hilbert spaces 
with a finite field, we must therefore drop the requirement for 
inner products and content ourselves with a plain vector space. 
Furthermore, in the absence of a notion of length, one must 
exclude zero vectors. 

Invertible Linear Maps: In actual QC, the dynamic evolu- 
tion of quantum states is described by unitary transformations 
which preserve inner products. As discrete quantum theory 
lacks inner products, the dynamic evolution of quantum states 
is described by any invertible linear transformation, i.e., by 
any linear transformation that is guaranteed never to produce 
the zero vector 

As an example, there are 16 linear (not necessarily in- 
vertible) functions in the space of 1-qubit functions. Out of 
these, six are permutations on the three 1-qubit vectors; the 
remaining all map one of the vectors to • which makes 
them non-invertible. Hence the space is quite impoverished 
compared to the full set of 1-qubit linear transformations in 
the Hilbert space. In particular, even some of the elementary 
unitary transformations such as the Hadamard transformation 
are not expressible in that space. Indeed the Hadamard matrix 
for the field of booleans maps the vector |+) to the zero vector 
and hence is not an acceptable 1-qubit transformation. 

Entanglement and Superdense Coding: Despite the re- 
striction to finite fields and the drastic reduction in the state 
space of qubits and their transformations, the theory built on 
the field of booleans has a definite quantum character. We 
present the superdense coding example from the paper of 
Schumacher and Westmoreland |[T| in Section lVIII-Cl and refer 
the reader to their paper for more details. 

V. Outline of Technical Development 

Our aims are to (i) develop a typed programming language 
that pins down the unique features of discrete quantum compu- 
tation and, (ii) to relate this language to that of conventional 
relational programming. Building on the work of Abramsky 
and Coecke (2], lH and of Selinger H, Q, both tasks can be 
achieved by distilling a logic and type system from a special 
category of sets and relations and by designing an appropriate 



syntax. That category (and hence the language) can be related 
to the conventional category of sets and relations to establish 
a connection with conventional relational programming and it 
can be connected to discrete quantum computation by showing 
that it is dagger compact closed. We proceed in the following 
two stages. 

Classical Computation: Many programming models of 
QC start with the A-calculus as the underlying classical 
language and add quantum features on top of it [20 1, I 2T1 . 
II22I . ||23]| . This strategy is natural given that the A-calculus 
is the canonical classical computational model. However this 
strategy complicates the development of quantum languages as 
it forces the languages to deal in fairly complicated ways with 
the implicit duplication and erasure of information in the clas- 
sical sublanguage. A simpler strategy that loses no generality 
is to build the quantum features on top of a reversible classical 
language: this keeps the language simple while still enabling 
the A-calculus constructs to encoded using two explicit op- 
erators for erasure and duplication if desired ll24l . Following 
this strategy we review the reversible language 11 that was 
recently developed by the first and third author and use it as a 
foundation for the quantum language in the remainder of the 
paper The language provides the symmetrical monoidal part 
of the required categorical structure. 

Monads and Arrows for Quantum Computation: We 
define a language TZIi by adding a layer on top of the classical 
reversible language 11 that provides sets and relations. The set 
construction can be expressed as a strong monad [25] over the 
language of classical observable values. However in order to 
express the required compact closure structure in the language, 
the Kleisli maps of this monad are abstracted in a first-order 
data type using an arrow fT4]. The language TUl is expressive 
enough to implement the superdense coding protocol and other 
quantum protocols. It differs from a conventional relational 
programming language in the semantics of the "union" oper- 
ation. 

VI. Classical Computation 

We introduce the language 11 for classical observable values. 
The language is defined over finite types: the only "computa- 
tions" in it are the isomorphisms between these types. In the 
context of Hilbert spaces, this language formalizes the notation 
used for the bases. In the categorical context, this language 
provides the underlying symmetric monoidal structure. 

A. Bases for Vector Spaces 

In the traditional presentation of quantum computing, a 
Hilbert space of dimension d is spanned by a set of d mutually 
orthogonal vectors. For example: 

« a quantum system Q2 composed of just 1 qubit is of 
dimension 2 and is spanned by {|0), |1)}; 

• a quantum system Q4 composed of 2 qubits is of dimen- 
sion 4 and is spanned by {|00), |01), |10), |11)}. 
The Hilbert space for the composition of two systems is 
spanned by basis vectors of the form \m ■ n) where m is a 
basis vector of the first system and n is a basis vector of the 



second system and the operation ( • ) is the composition of 
the two labels. For example: 

• the composition of Q2 and (written Q2 ® Qa) is of 
dimension 8 and is spanned by: 

{|0-00), 10-01), 10-10), 10-11), 11-00), 11-01), 11-10), 11-11)} 

• the composition of (54 and Q2 (written ® Q2) is also 
of dimension 8 and is spanned by: 

{|00-0), 100-1), 101-0), 101-1), 110-0), 110-1), Ill-O), lll-l)} 

Clearly the compositions Q2®Qi and Qi®Q2 are isomorphic 
and the distinction between them is typically blurred. 

B. Syntax and Type System 

Instead of using bits as the labels for the bases, we use 
structured finite types built using the empty type, the unit type, 
the sum type, and the product type. Furthermore instead of 
silently identifying systems like Q2 ® Qa and Q4 (g) Q2, we 
include an explicit set of operators between these finite types to 
witness the isomorphisms. Specifically, we have the following 
language of classical observable types b and values v: 

b ::= 0|l|5 + 5|6x6| bool 

V ::= \ left V \ right v \ {v,v) \ 1 \ F 

The types are finite as they will be used to define the 
dimensions of our vector spaces. The type is the empty 
type containing no inhabitants. The type 1 has exactly one 
inhabitant called (). The type 61 + 62 is the disjoint union of 61 
and 62 whose elements are appropriately tagged values from 
either type. The type 61 * 62 is the type of ordered pairs whose 
elements are coming from 61 and 62 respectively. Although the 
type of booleans is expressible as 1 + 1 we add it as a primitive 
type because it corresponds to the type of scalar field values 
which play an important role in the language. The type system 
is summarized below: 

h u : 61 h w : 62 

h : 1 h left v:bi+b2 h right v : 61 + 62 
h Hi : 61 h t>2 : 62 

h {vi,V2) : 61 X 62 h T : 600/ h F : 600/ 

The following isomorphisms are sound and complete for the 
given finite types [261 : 
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X63) 
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X 62) X 63 
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) X63 




ih 


X 63) + (62 
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1 



X63) 



The forward c vi V2 transitions are given below: 



bi o 62 ^2 ^ &3 
h ^ 63 

bi -H- 63 b2 bi bi -H- 63 62 bi 
(61 + 62) <H- (63 + 64) (61 X 62) <H- (63 X bi) 

We introduce primitive operators corresponding to the left- 
to-right and right-to-left reading of each isomorphism. We 
gather these operators into the table below. 



id 
>+ 
x+ 
^+ 

Xx 

^0 
-< 



b^b 

+ b^b 

61 + 62 ^ 62 + ^1 
61 + (62 + 63) ^ (?>i + b2) + 63 

1 X 6 6 

&i X 62 &2 X 61 
61 X (62 X 63) (61 X 62) X 63 

X 6 o 

(61 + 62) X 63 (61 X 63) + (62 X 63) 

bool -H- H- 1 



id 
<+ 
x + 

^+ 

=^x 
Xx 

^x 
^0 
>- 

5 



Each line of this table is to be read as the definition of two 
operators. For example corresponding to the identity of x 
isomorphism we have the two operators ^ x : 1 x b -H- b and 
<x- b ^ 1 X b. 

Now that we have primitive operators we need some means 
of composing them. We construct the composition combinators 
out of the closure conditions for isomorphisms. Thus we have 
one sequential composition combinator o and two parallel 
composition combinators, one for sums © and one for pairs ®. 



ci : bi O 62 C2 : 62 ^ ^3 
(ci o C2) : 61 -H- 63 



Ci : &i O 63 C2 : 62 O bj 
(ci (g) C2) : (61 X 62) <^ {ba x 64) 

Definition 6.1: (Syntax of 11) We collect our combinators, 
types and values to get the definition of our language for 
isomorphisms, which we will refer to as 11: 
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::= 0\l\b + b\ bxb 


bool 
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■■■■= 1 left V 


right v 


{v,v) 


F 1 T 
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■■■■= x+ 1 ^+ 


^+ 1 


>+ 1 


< + 




1 Xx 1 ^x 


1 ^x 1 


^x 1 


^x 




1 ^0 1 ^0 1 


K 1 >- 1 


id € 


1 5 


c 


::= iso 1 c c 


c c 1 c 


©c 





Given a program c : 61 -H- 62 in 11, we can run it either in 
the forward direction by applying it to a value vi : bi or run 
it in the opposite direction by applying it to a value V2 ■ 62. 
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Jpff 77i 
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( Ipff ( rinhf Vn)^ 
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right () 




(left 0) 
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{right 0) 




F 



ci : 61 O 63 C2 : &2 ^ bj 
(ci ©C2) : (61-^-62) ^ (63-1-64) 



Since there are no values that have the type 0, the reductions 
for the combinators )>=+, =^+, ^0 and ^0 omit the impossible 
cases. The semantics of the other combinators is straightfor- 
ward. Composition combinators are defined as follows: 

Ci Vi !->■ V2 

(ci © C2) {left vi) !->• left V2 

C2 Vi ^ V2 



(ci © C2) {right ui) 1— > right V2 

Ci Vi H> '(;3 C2 1)2 I— f 4 Ci Vi i-^ V C2 V ^ V2 



(ci (g) C2) (t;i,w2) 1-^ (^^3,^4) (ci o C2) W2 

For the inverse direction, it is a simple matter to estabhsh 
the following property. 

Proposition 6.2: For every 11 program c, such that c v ^ 
v', we can construct its adjoint such that c^v' i-> v. 

Proof: We can construct the required c' by replacing 
every primitive isomorphism with its dual. For sequential 
composition we have (ci 002)^ = ocj and for parallel com- 
position we have (ci(8)C2)^ = c|(8)C2 and (ci®C2)^ = c|®C2. 



C. Dagger Symmetrical Monoidal Structure 

The language 11 is rich enough to express the required 
dagger synmietric monoidal (but not the compact close struc- 
ture) structure. The proof for this is straightforward and we 
content ourselves with a brief outhne. The language 11 can be 
interpreted as category whose objects are the base types 6 and 
whose morphisms the combinators c. Identity morphisms are 
given by id and associativity follows from the composition of 



isomorphisms, thus establishing the required properties for a 
category. 

Further 11 is a monoidal category {b, x, 1) where x is the 
tensor operation, 1 is the identity object and the required 



and 



natural transforms a, A and p are given by 
X X o ;>= X respectively. The required "pentagon" and "triangle" 
axioms follow from the definitions of the these isomorphisms. 

Braiding is provided by x ^ at the appropriate types and it 
satisfies the "hexagon" axioms. Further that x x at fei x 62 is 
the same as Xx^ at 62 x 61 establishing symmetry. 

Finally, proposition 16.21 tells us that for every morphism c 
there exists its adjoint establishing that 11 is a dagger 
symmetric monoidal category. 

VII. The Language Ten 

The core language for discrete quantum computing is ob- 
tained by adding a layer on top of the reversible core presented 
in the previous section. The additional layer is that of vectors 
and linear maps or equivalently that of sets and relations. 

Definition 7.1: (Syntax of Core 7^,11) 



b 

V 

s 
r 



= ...\Sb\bRb 
= ... |. 

= I {w} I s ttl s 
— arr iso \ r ^ r | second r 
strength \ state s | ?7f, | e 



The language TZIV extends 11 as follows. The set of types 
is extended with the type S* 6 of sets of values and the type 
&i R 62 of relations between sets of values of type bi and sets 
of values of type 62- The set of values is extended with sets s 
which can either be the empty set 0, a singleton set {u}, or 
the exclusive union si tbi S2- The exclusive union of sets si 
and S2 is the union of all the elements that are not in common. 
This union reflects the cyclic nature of the underlying finite 
field as was explained in Section |IV] In more detail, we saw 
for example that adding the vector f|0) + t|1) to the vector 
t|0) + t|1) produces the vector t|0) + f|1) where the 
two occurrences of the component |1) canceled each other 
Expressed as sets, we would have that {1} l±l {0, 1} is equal 
to {0}. 

The layer of sets corresponds to a strong monad with the 
singleton set as the unit and the folding of the exclusive 
union as the bind operation. The Kleisli arrows of this monad 
are functions that map values to sets. In order to express 
these functions themselves as relations, we provide explicit 
arrow operators to construct them. The language therefore 
includes a separate syntactic category of relations which is 
used to build the required Kleisli maps using arrow primitives. 
The operator arr lifts any isomorphism from the underlying 
reversible language to a relation on sets. The operator ^ 
sequences two relations and the operator second applies a 
relation to the second component of a pair leaving the first 
component unchanged. These three operators are the minimal 
core operators for any arrow language. In addition, we must 
include enough structure to provide a compact closed category: 
the operator strength is used to express the tensor product of 



sets; the operator state identifies a state of type S b with an 
arrow (relation) of type 1 R b f3l; finally the pair of opera- 
tors rib and e are the unit elimination and introduction that are 
required for any compact closed category. The operation 77 is 
indexed by a type b. 

We present the type system below. The remainder of this 
section explains the semantics of core TZIi formally and 
illustrates its use in several examples: 

\- V -.b h si : S b \- S2 ■ S b 

h : 5* 6 h{v}:Sb h si W S2 : 5* & 

\- iso -.b ^ b \- ri : bi R b2 h r2 : 62 63 



h arr iso : b R b h ri r2 '■ bi R 63 

\-r:biRb2 h s : S b 



I- second r : (b x bi) R {b x 62) ^ state s : I R b 
h strength : (bi x S ^2) R {bi x 62) 



h Tyfc : 1 i? (6 X 6) h e : (6 x &) i? 1 
A. Semantics 

A program in TZIl consists of the application of a relation r 
to a set s producing a resulting set s' . We model this evaluation 
using two applications: the application r @ s applies the rela- 
tion r to the set s and the application r@v applies the relation 
r to an individual set element v. The first application @ 
simply iterates down the structure of the set until it finds 
individual elements which can be processed using @ : 

r"¥0 ^ 

r'@{v} M- r@v 

r^(si 1+1 S2) ^ (r^si) 1+1 {r'@S2) 

The application @ applies a relation to a single value v and 
returns the set of possible values that are related to v: 

[arr iso)@v 1— >■ {iso{v)} 

{ri r2)@v r2'@{ri@v) 

{second r) @ {vi,V2) 1— )> strength @ {vi,r @ V2) 

strength @ (u, 0) 

strength @{vi,{v2}) {(wi,W2)} 

strength @ (v, s I 1+) S2) H> strength @ {v, si) W 

strength @ (w, S2) 

{state s)@{) H- s 

Vb@{) H- 1+J,; u,)}: whereu^ e 6 

s@{v,v) H- {()} 

e@{v,v') ^ 9 if V ^v' 

When applied to a value v, the relation arr iso simply 
applies the underlying isomorphism to the value which is 
wrapped in a singleton set. Applying the sequence ri and r2 
to a value v applies ri to v to produce a set which is passed 
to r2- To apply a relation r to the second component of a 
pair {vi,V2), we first produce {vi,r@V2) and then use the 
tensorial strength of the monad to push the pair construction 
through the resulting set as shown in the three rules for 
strength. The relation state s maps the unit value () to the 
given set s. When applied to the unit value (), the relation r/i, 



produces a superposition of all pairs of values Vi where Vi 
is an element of the type b (which is finite). The relation e 
maps a pair of equal values to the singleton set containing () 
and maps a pair of different values to the empty set. As the 
semantics for r/b and e shows, these constructs exploit the 
fact that the underlying language of classical values is based 
on finite values (and hence values that can be compared for 
equality and enumerated). In order to accommodate sets whose 
elements are themselves sets or relations, it is necessary to 
have a first-order representation of relations using the arrow 
constructors instead of using the higher-order Kleisli maps of 
the monad. 

It is a routine task to verify that the evaluation rules preserve 
the types. 

Proposition 7.2: If h r : 61 i? 62 and ^ s : S hi then 
h r~@s : S b2. 

B. Derived Relations 

The language 7^11 is surprisingly expressive: using standard 
categorical constructions, it can express higher-order func- 
tions, currying, uncurrying, recursion, adjoints, dot products, 
and outer products. We illustrate these derived relations. The 
superdense coding protocol in Section IVIII-CI uses the 
construction below. 

We start by defining for every relation r a relation first r 
that applies r to the first component of a pair: 

first : a R b ^ {a X c) R (b X c) 
first r = arrX y 3^ second r ^> arrX y 

Next we present currying and uncurrying and use them to 
turn any set of pairs to a relation: 

curry : ((c x a) R b) ^ {c R {a x b)) 
curry r = arr °Xx) second rja ^ 
arr ^ x 3S> first r arr x x 

uncurry : {c R {a x 6)) — > ((c x a) R b) 
uncurry r — first r 3S> arr ( x x o ^ x ) 

S2r : S (a, b) (a R b) 
S2r s = arr =^x uncurry (state s) 

We can also define a trace operation to model recursion 
from cyclic sharing 117|: 

trace : {{a x c) R (b x c)) (a R b) 
trace r — arr =<; x Sl> first rjc ^ 

arr (x^o ^x) ^ first r :s> arr ^x 
^ second e arr ( x x o ^ x ) 

Finally, we can also define an adjoint for every relation 
and use it to define a costate relation that matches a given 
state. We can combine a state and a costate in two different 
ways to simulate the dot product and the outer product 



constructions lH]. 

adjoint : {a R b) (b R a) 
adjoint r — arr ^ first rja ^ 

first (second r) arr ^x 3S> 
second e ^ arr (XyO )>= x ) 

costate : 5 a — > (a i? 1) 
costate s — adjoint (state s) 

(. I .) : Sa^Sa^(lRl) 
(si I S2) = state si ;:g> costate S2 

|.)(.| : Sa^Sa^(aRa) 
■5i)('*2| — costate S2 state si 

The result of the dot product is a relation of type 1 R 1 which 
corresponds to a scalar ||3|. 

C. Interpretations 

There are two isomorphic ways of thinking about TZIl: a 
value of type S b could be viewed as a set of values of 
type 6 or as a vector which maps values of type b to scalars 
in the field of booleans. For example, s : S (hool x hool) — 
{(f, F)}ttl {(T, t)} denotes the set {(f, f), (t, t)} in the first 
interpretation and denotes the vector: 

(f,f) /t\ 
(f,t) f 
(t,f) f 

(t, t) yi ) 

in the second interpretation. In the vector notation, the labels 
on the left enumerate all the values of the given type. A value 
is present if the corresponding entry in the vector is T and 
absent if the corresponding entry is F. To avoid clutter we 
will assume a fixed preferred ordering of the labels on the left 
and omit them. 

Similarly a value of type &i R 62 can be viewed as a relation 
mapping sets of values of type &i to sets of values of type 62 
or as a linear map which given our preferred ordering could 
be represented as a matrix. For example, the values: 

ri,r2 '■ hool R hool 

n = S2r ({(F,F)}W {(T,F)}a {(T,T)}) 
r2 = S2r ({(F,F)} W {(F,T)} a {(T,F)}) 

denote the relations {(f, f), (t, F), (t, t)} and 
{(f, f), (f, t), (t, f)} in the first interpretation, and the 
matrices: 

r.= (^ r2=(^ ^\ 

in the second interpretation. In the matrix notation, the 
columns are implicitly indexed with F and T from left to right, 
and the rows are implicitly indexed with F and T from top to 
bottom. An entry is T if the pair (column-label, row-label) is 
in the relation and F otherwise. 



In a conventional setting, the composition of ri and 7'2 
would be: 

r2 = {(F, F), (F, T), (T, F), (T, T)} 

This is the semantics one gets in conventional relational 
programming, i.e., pure Prolog and backtracking monads |27|, 
||28 1 which indeed can be implemented and formalized using 
sequences or sets 1291 . 

However in 7^11, the composition of these two same rela- 
tions is {(f, f), (f, t), (t, t)} because the pair (T, F) can be 
produced in two different ways which cancel due to interfer- 
ence. It is perhaps more intuitive to compute the composition 
in the world of matrices: 



TVF TVT \ _ / T F 
TVF TVF / ^ V T T 



We can give similar interpretations for each of the TZIl 
constructs. In the world of matrices, second r denotes the 
matrix produced by the tensor product of the identity matrix 
and the matrix corresponding to r. The construct state is a 
no-op: it allows us to view a vector of size n as a matrix of 
dimensions nxl. We illustrate the matrix denoted by strength 
at the type (bool x S hool) R {boot x bool) 

T F T F F F f\ 

F T T F F F F 

F F F F T F T 

F F F F F T T/ 



F 
F 



(T 



F F 



T) 



The rows are indexed from top to bottom by (f.f), (f,T), 
(t, f), and (t, t). The columns are indexed from left to right 
by (F, 0), (F, {F}), (F, {T}), (F, {F} W {t}), (T, 0), (T, {f}), 
(T,{T}), and (T,{F} W {T}). 

The matrices corresponding to and e are column and row 
vectors with entries T at the diagonal elements. For example, 
rjhooi ■ 1 R {bool X bool) and e : {bool x bool) R 1 are: 

fT\ 

F 
F 
\jj 

As a special case, we note that relations of type 1 i? 1 denote 
matrices of dimension 1x1, i.e., scalars. We interpret the 
matrix (t) as the scalar T and the matrix (f) as the scalar F. 

Lemma 7.3: The semantics of 7^11 is sound with respect to 
both the relational interpretation and the vector interpretation. 

D. Dagger Compact Closed Structure 

We can now verify that the arrow fragment of the language 
7?.n forms a dagger compact closed category. The objects 
of the category are the types S b and the morphisms are 
the relations bi R b2- These form a category since we have 
identity morphisms due to arr id and the composition of 
relations is associative due to the associativity of the matrix 
multiplication in the model (soundness). The category has the 
required monoidal structure with S {bi x 62) being the product 
for the objects S bi and S 62 and 5 1 is the identity object 



for tensors and the required natural transforms a, A and p are 
arr lifted forms of their 11 equivalents. Similarly braiding and 
symmetry are preserved by lifting x x from 11. 

The dagger structure on relations is given by the contravari- 
ant functor adjoint : 5i i? 62 — >■ ^2 which we constructed 
in Section IVII-BI and is easily verified to be involutive. 

Finally the required autonomous structure is given by choos- 
ing the dual of S" 6 to be the object S b itself. The required 
unit and counit are given by 77 and e in 7^11. The required 
adjunction triangles become the same as the right and left 
duals on objects coincide. Thus the arrow fragment of 7^11 
forms a dagger compact closed category. 

VIII. Using VM for Quantum Computation 

Two more things are needed to use TUl for quantum 
computation. First, the zero vector (which is denoted by 0) 
should never be encountered during computation as it does 
not correspond to a valid quantum state. Second, we must 
define a notion of measurement. 

A. Evolution of Quantum States 

As Schumacher and Westmoreland fV\ explain, the evolution 
of quantum states in DQC2 only requires that the linear 
maps be invertible. We have already established that every 
linear map has an adjoint but these adjoints do not necessarily 
coincide with the inverses. Consider the matrix: 



T T 
T T 



The adjoint of this matrix is itself but when multiplied by 
itself in the field of booleans it produces the zero matrix. In 
other words, any TUl expression that can denote this matrix 
should not be allowed. All the expressions in the syntactic 
category r in TZIi denote invertible matrices except for e. 
This suggests a possibility for tracking the uses of e using an 
extended type system to syntactically guarantee that a given 
expression denotes an invertible transformation. 

B. Measurement 

In standard QM, the postulate of measurement requires the 
notion of an orthogonal projection. However, as explained in 
Section IrVl discrete quantum theories lack inner products and 
hence lack a standard notion of orthogonality. Nevertheless 
it is possible to define a sensible notion of measurement as 
explained by Schumacher and Westmoreland fll. The key 
technical observations are (i) any set of linearly independent 
vectors can form a basis, and (ii) it is possible to define dual 
vectors for a given basis. 

In more detail, assuming the vector interpretation of 7^11, 
a measurement (or an observable) corresponds to a basis 
{.si, S2, . . .} where the collection of values Si denote linearly 
independent vectors. As discussed in Section IIVI the space 
of 1-qubit vectors contains just three vectors with any pair 
forming a basis. For each choice of basis, we associate an 



observable: 



X-basis = {{t},{f}W {t}} 
Y-basis = {{f} l±) {t},{f}} 
Z-basis = {{f},{t}} 

We then associate a basis dependent dual Ji to each vector s; 
such that (s7 | Sj) denotes the scalar T if and only if i — j. 
For example, the dual vectors for the vectors in the X-basis 
are: 

X-dualx = {{f}1±) {t},{f}} 

The simplest way then to extend TZIi is to add one construct 
measure that takes two arguments: a vector and a set of dual 
basis vectors along which the vector should be measured. 
For example, to measure the vector {t} in the X-basis, we 
would write: measure {t}{{f} l+l {t},{f}} The result of 
measurement is any dual basis vector that can possibly match 
the given vector selected at random. In the above example, 
only the first dual basis vector can possibly match the given 
vector which means that the result is deterministic. 

A remarkable feature of the categorical approach to QM is 
that the semantics of measurement can be expressed as follows 
in 7^^: 



where dualhasis = { 



measure s 



{di,d2, .-•} = di if {di I s) denotes (t) 



As explained in Sections IVII-BI and IVII-CI the dot product is 
expressible in TZIi and the relations of type 1 i? 1 are identified 
with scalars. If more than one dual vector produces the scalar 
T then the result of measurement is non-deterministic with the 
dual vector di picked at random. 

It should also be possible to exploit the fact that our 
categories have biproducts to extend the language in a richer 
way by allowing measurements to occur in the middle of 
computation [3J. 

C. Example: Superdense Coding 

The superdense coding example presented by Schumacher 
and Westmoreland |r| can be directly implemented in TZIi. 
Alice and Bob initially share an entangled state, represented by 
{(f, f)} I+I {(t, t)}. Depending on Alice's choice of sending 
numbers to 3, Alice applies one of the following operations 
on the first bit: 

a/ice = arr id 

alicei — arr o x x o 5) 

alice2 = S2r {state ({(f,f)}W {(T,F)}a {(t,t)})) 
alices = alicei ^ alice2 

By measuring in a particular dual basis described below. 
Bob will deterministically set a different dual vector for each 
possible operation that Alice could have performed. The entire 
example is then written as follows: 

measure {{first a/ice„) W ({(f, F)} 1+) {(t,t)})) dualhasis 



{(F,T)}W {(T,F)}W {(T,T)}, 

{(F,F)}W {(T,F)}W {(T,T)}, 

{(F,T)}W {(T,F)}, 

{(F,F)}W {(T,T)}} 

Interestingly, it is possible to write this same example in 
conventional relational programming. The fact that intermedi- 
ate results are accumulated using the standard union instead of 
the exclusive union is apparent in the results. However because 
this computation does not use significant intermediate steps, 
the only values that need to be treated specially are the final 
results. Indeed if the results that appear an even number of 
times are removed, then the Prolog execution performs the 
superdense coding exactly. 

The complete program and its execution are below: 



r ( false , false ) . 
r (true , true ) . 

s ( false , true) . 
s (true , false) . 

u ( false , false) . 
u ( false , true) . 
u ( true , true) . 

V ( false , false) . 
V ( false , true) . 
v(true, false) . 

rd( false , true ) . 
rd(true , false) . 
rd(true , true) . 

sd( false , false) 
sd(true , false) . 
sd(true , true) . 

ud( false , true ) . 
ud (true , false ) . 

sdcoding(N,M) :- 



vd( false , false ) 
vd(true , true) . 

eq( false , false ) 
eq (true , true) . 

id ( false , false ) 
id (true , true) . 

g ( false , true) . 
g ( true , false) . 

k ( false , false ) . 
k ( true , false) . 
k ( true , true ) . 



gk(X,Y) 



g(x,z) ,k(z,Y) . 



alice(0 ,x, Y) : — 
alice(1 ,x, Y) : — 
alice(2 ,x, Y) :- 
alice(3 ,x, Y) : — 



id(x,Y) . 
g(x,Y) . 
k(X,Y) . 

gk(x,Y) . 



r(X,Y) ,alice(N,X,B) .measure (( B , Y) ,M) . 



measure ( ( SI , S2 ) ,0) 
measure ( ( SI , S2 ) ,1 ) 
measure ( ( SI , S2 ) ,2) 
measure ( ( SI , S2 ) ,3) 



- rd(Bl,B2) ,dotP((Bl,B2) ,(S1,S2) ) . 

- sd(Bl,B2) ,dotP((Bl,B2) ,(S1,S2) ) . 

- ud(Bl,B2) ,dotP((Bl,B2) ,(S1,S2) ) . 

- vd(Bl,B2) ,dotP((Bl,B2) ,(S1,S2) ) . 



dotP((Bl,B2) ,(S1,S2) ) 



28 ?- sdcoding (0, X) . 
X = 1 
X = 3 
X = 
X = 1 
X = 3 . 

29 ?- sdcoding (1, X) . 
X = 
X = 1 
X = 2 
X = 
X = 2 . 

30 ?- sdcoding (2, X) . 
X = 1 ; 



eq(Bl , SI) ,eq(B2 , S2 ) . 



X = 3 

X = 

X = 2 

X = 

X = 1 

X = 3 . 



31 ?- sdcoding(3,X) 
X = 1 
X = 3 
X = 
X = 1 
X = 2 
X = 
X = 2 . 



IX. Conclusion 

We have distilled a programming model for the discrete 
quantum theory over the field of booleans recently introduced 



by Schumacher and Westmoreland 0]. The model is expressed 
in a small calculus 7^11 with formal type rules and semantics. 
The language TZIl is directly inspired by the computational 
structures of quantum mechanics previously identified by 
Abramsky and Coecke and Selinger. The semantics of TZIl is 
sound with respect to both a relational model with exclusive 
unions or a vector space model over the field of booleans. 
Computationally 7^11 is a relational programming language 
surprisingly similar to traditional relational programming lan- 
gauges like Prolog with the significant difference that where 
Prolog would accumulate all possible answers with repetitions, 
TZIl allows interference between possible answers. 

A natural goal of our research is to unravel the mathematical 
underpinnings behind the power of quantum computation. By 
analyzing models such as those presented in this paper, we 
hope we will be able to shed some light on this fundamental 
issue. Deutsch's quantum algorithm, which establishes in a 
single measurement whether a boolean function is constant 
or balanced, makes use of interference and quantum paral- 
lelism. Although discrete quantum computing over the field 
of booleans, DQC2, contains both these properties, Deutsch's 
algorithm cannot be efficiently implemented within this model. 
This example illustrates the non-trivial character of the quest to 
disentangle the power of quantum computation. In more detail, 
having eUminated much of the structure of actual quantum 
mechanics, the connection between DQC2 and conventional 
relational programming singles out the use of exclusive union 
as the source of all the "quantum-ness". Assuming that the 
exclusive union would be performed in "constant time" by a 
quantum computer, more work is needed to establish whether 
it is possible to write more efficient algorithms using the 
exclusive union. 
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